Ab initio description of the exotic unbound 7 He nucleus 
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The neutron rich exotic unbound 7 He nucleus has been the subject of many experimental investiga- 
tions. While the ground-state 3/2~ resonance is well established, there is a controversy concerning 
the excited l/2~ resonance reported in some experiments as low-lying and narrow (Er~1 MeV, 
r< 1 MeV) while in others as very broad and located at a higher energy. This issue cannot be 
addressed by ab initio theoretical calculations based on traditional bound-state methods. We in- 
troduce a new unified approach to nuclear bound and continuum states based on the coupling of 
the no-core shell model, a bound-state technique, with the no-core shell model/resonating group 
method, a nuclear scattering technique. Our calculations describe the ground-state resonance in 
agreement with experiment and, at the same time, predict a broad l/2~ resonance above 2 MeV. 



PACS numbers: 21.60.Dc,24.10.Cn,25.10.+s,27.20.+n 



The 7 He nucleus is an exotic system of three 
neutrons outside a 4 He core with a particle- 
unstable JT=3/2 _ 3/2 ground state (g.s.) lying 
at 0.430(3) MeV [H, above the threshold of a neutron 
and 6 He, which in turn is an exotic Borromean halo 
nucleus. While there is a general consensus on the 5/2~ 
resonance centered at 3.35 MeV, which mainly decays to 
a+3n [H, discussions are still open for the other excited 
states. In particular, the existence of a low-lying 1/2" 
state at about 1 MeV has been advocated by many ex- 
periments 4-8] (most of them using knockout reactions 
with a 8 He beam on a carbon target), while it was not 
confirmed in several others [9 14]. This contradictory 
situation arises from the main experimental difficulty 
of measuring the properties of excited states in 7 He 
in the presence of a three-body background, coming 
from the particle decay of 7 He and from the outgoing 
particle involved in the reaction used to produce 7 He. 
The presence of a low-lying 1/2" state has also been 
excluded by a study on the isobaric analog states of 7 He 
According to this latter work, a broad l/2~ 
at -3.5 MeV with a width T-lOMcV fits 
the data the best. Neutron pick-up and proton-removal 
reactions 11, [l2| suggest instead a l/2~ resonance at 
about 3 MeV with a width T«2 MeV. 

The l/2~ resonance controversy cannot be addressed 
by ab initio theoretical calculations based on traditional 
bound-state methods such as the Green's function Monte 
Carlo (GFMC) 0, the no-core shell modeWNCSM) [T3] 
or the Coupled Cluster Method (CCM) [3S]. The 
complex CCM was recently applied to He isotopes, but 
only the g.s. of 7 He was investigated [211 ] . 

In this Letter, we address the low-lying resonances of 
7 He within the no-core shell model with continuum (NC- 
SMC), a new unified approach to nuclear bound and 
continuum states based on the coupling of the NCSM 
with the no-core shell mo del/ resonating group method 
(NCSM/RGM) (22M27I]. In this approach, we augment 
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the NCSM/RGM ansatz for the A-body wave func- 
tion [23| by means of an expansion over A-body NCSM 
eigenstates \AXJ^T) according to: 
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where the (A — a, a) binary-cluster channel channel states 



\<r T ) = [(\A~a aJ^T^a a 2 I?T 2 )f T ^ 

(J'T) S (r-r A -a,a) 



X Yi(fA-a,a) 



rr A - 



(2) 



are labeled by the collection of quantum numbers v — 
{A — a ail^T^a ail^T^sl} and rA-a,a is the inter- 
cluster relative vector. The NCSM sector of the basis 
provides an effective description of the short- to medium- 
range A-body structure, while the NCSM/RGM clus- 
ter states make the theory able to handle the scattering 
physics of the system. The discrete, c\, and the contin- 
uous, 7i/(r) unknowns of the NCSMC wave functions are 
obtained as solutions of the following coupled equations, 
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where Xv( r ) are the relative wave functions in the 
NCSM/RGM sector when working with the orthogo- 



nalized cluster channel states [23J. The NCSM sector 
Hncsm of the Hamiltonian kernel is a diagonal matrix of 
the NCSM energy eigenvalues, while rl is the orthogonal- 
ized NCSM/RGM kernel [23j]. The coupling between the 
two sectors is described by the overlap, g\ v { r ), and hamil- 
tonian, h\ v (r), form factors respectively proportional to 
the (AXJ^TIA^^) and (AAJ 7r T|ffX|*^ T ) matrix 
elements. We solve the NCSMC equations by applying 
the coupled- channel microscopic R-matrix method on a 
Lagrange mesh jjij ]. Further details on the formalism 
will be given elsewhere [29[ . 
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Eg. s . [MeV] 


An 

4 He 


b He 


7-TT 

7 He 


NCSM A max =12 


-28.05 


-28.63 


-27.33 


NCSM extrap. 


-28.22(1) 


-29.25(15) 


-28.27(25) 


Expt. 


-28.30 


-29.27 


-28.84 



TABLE I: Ground-state energies of 4 ' 6 ' 7 He in MeV. An expo- 
nential fit was employed for the extrapolations. 



We begin by presenting NCSM calculations for 6 He 
and 7 He that will serve as input for the subsequent 
NCSM/RGM and NCSMC investigations of 7 He. In 
this work, we use the similarity-rnormalization-group 
(SRG]evolved chiral N 3 LO NN potential of 

Refs. [34l 35 1 . For the time being, we omit both induced 
and chiral initial three-nucleon forces, and our results de- 
pend on the low-momentum SRG parameter A. However, 
for A = 2.02 fm , we obtain realistic binding energies 
for the lightest nuclei, e.g., 4 He and, especially important 
for the present investigation, 6 He (see Table [TJ) . Conse- 
quently, this choice of TV TV potential allows us to perform 
qualitatively and quantitatively meaningful calculations 
for 7 He that can be compared to experiment. Except 
where differently stated, all results shown in this work 
have been obtained with an harmonic oscillator (HO) 
N max =12 basis size and frequency MeV. 

The variational NCSM calculations converge rapidly 
and can be easily extrapolated. At N mlix —12 (our 6 ' 7 He 
limit for technical reasons), the dependence of the 6 He 
g.s. energy on the HO frequency is flat in the range 
of M7 ~ 16—19 MeV. In general, when working within 
an HO basis, lower frequencies are better suited for the 
description of unbound systems. Therefore, we choose 
h£l=16 MeV for our subsequent calculations. Extrap- 
olated g.s. energies with their error estimates and the 
TV max =12 results are given in Table HI Calculated 6 He ex- 
citation energies for basis sizes up to A^ max =12 are shown 
in Fig.[TJ The 6 He is weakly bound with all excited states 
unbound. Except for the lowest 2 + state, all 6 He excited 
states are either broad resonances or states in the con- 
tinuum. We observe a good stability of the 2f state 
with respect to the basis size of our NCSM calculations. 
The higher excited states, however, drop in energy with 
increasing iV max with the most dramatic example being 
the multi-foil 0^ state. This spells a potential difficulty 
for a NCSM/RGM calculations of 7 He within a n+ 6 He 
cluster basis as, with increasing density of 6 He states at 
low energies, a truncation to just a few lowest eigenstates 
becomes questionable. 

For the 7 He, the NCSM predicts the g.s. unbound in 
agreement with experiment. However, the resonance en- 
ergy with respect to the 6 He+n threshold appears over- 
estimated. Obviously, it is not clear that the ad hoc 
exponential extrapolation is valid for unbound states. In 
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FIG. 1: (color online). Dependence of 6 He excitation energies 
on the size of the basis N max . 



addition, no information on the width of the resonance 
can be obtained from the NCSM calculation. We can, 
however, study the structure of the 7 He NCSM eigen- 
states by calculating their overlaps with 6 He+n cluster 
states, which are related to g\ v (see Eq. j2])), and the 
corresponding spectroscopic factors summarized in Ta- 
ble [Hi Overall, we find a very good agreement with the 
VMC/GFMC results as well as with the latest experi- 
mental value for the g.s. 0. Interesting features to no- 
tice is the about equal spread of l/2~ between the + and 
2^~ states. We stress that in our present calculations, the 
overlap functions and spectroscopic factors are not the 
final products to be compared to experiment but, on the 
contrary, inputs to more sophisticated NCSMC calcula- 
tions. 



7 He r 


6 He- 


-nilj) 


NCSM 


CK 


VMC 


GFMC 


Exp. 


3/2- 


0+ 


"Pf 


0.56 


0.59 


0.53 


0.565 


0.512(18) [2] 
0.64(9) [36] 
0.37(7) [11] 


3/2- 


2f 


-v\ 


0.001 


0.06 


0.006 






3/2- 


2f 




1.97 


1.15 


2.02 






3/2- 


2+ 

Z 2 


-v\ 


0.12 




0.09 






3/27 


2+ 

Z 2 


-p§ 


0.42 




0.30 






1/2- 


0+ 




0.94 


0.69 


0.91 






1/2- 


2+ 


-Pi 


0.34 


0.60 


0.26 






1/2- 


2+ 

z 2 


-p§ 


0.93 











TABLE II: NCSM spectroscopic factors com par ed to Cohen- 
Kurath (CK) [37| and VMC/GFMC HEIHlai] calculations 
and experiment. The CK values should be still multiplied by 
A/(A—1) to correct for the center of mass motion. 
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FIG. 2: (color online). Dependence of the NCSM/RGM (blue 
lines) and NCSMC (red lines) 6 He + n diagonal phase shifts 
of the 7 He 3/2~ g.s. on the number of 6 He states included in 
the binary-cluster basis. The short-dashed, dashed, and solid 
curves correspond to calculations with the 6 He + g.s. only, 
+ ,2 + states, and + ,2 + ,2 + states, respectively. 



The NCSM/RGM calculations for the n+ 6 He system 
presented in the following were obtained by including 
up to the three lowest eigenstates of 6 He, i.e., 0+,2+, 
and 2 J. These results will be compared to NCSMC 
calculations, which couple the above n+ 6 He binary- 
cluster states to the 6 lowest negative parity NCSM 
eigenstates of 7 He (3/27/ , 1/2", 5/2~, 3/2 J, 3/2,7, 3/2^ ) 
as well as the four lowest 7 He positive-parity eigenstates 
(1/2+5/2+ 3/2+, 5/2+). 

First, in Fig. [21 we study the dependence of the 3/2~ 
g.s. diagonal phase shifts on the number of 6 He eigen- 
states included in the NCSM/RGM (blue lines) and NC- 
SMC (red lines) calculations. The NCSM/RGM calcu- 
lation with the 6 He target restricted to its g.s. does not 
produce a 7 He 3/2~ resonance (the phase shift does not 
reach 90 degrees). A 2 P 3 / 2 resonance does appear once 
the 2+ state of 6 He is coupled, and the resonance posi- 
tion further moves to lower energy with the inclusion of 
the second 2+ state of 6 He. On the contrary, the 2 P%/2 
resonance is already present in the NCSMC calculation 
with only the g.s. of 6 He. In fact, this NCSMC model 
space is already enough to obtain the 7 He 3/2~ g.s. res- 
onance at about 1 MeV above threshold, which is lower 
than the NCSM/RGM prediction of 1.39 MeV when three 
6 He states are included. Adding the 2+ state of 6 He gen- 
erates a modest shift of the resonance to a still lower 
energy while the second 2+ state of 6 He has no signif- 
icant influence (Fig. [2J panel (b)). We further observe 
that the resonance position in the NCSMC calculation is 
lower than the NCSM/RGM one by about 0.7 MeV. This 
difference is due to the additional correlations brought 
by the 7 He eigenstates that are coupled to the n+ 6 He 
binary-cluster states in the NCSMC and that compen- 
sate for higher excited states of the 6 He target omitted in 




FIG. 3: (color online). NCSM/RGM (a) and NCSMC (b) 
6 He + n diagonal phase shifts (except 6 P%/2, which are eigen- 
phase shifts) as a function of the kinetic energy in the center of 
mass. The dashed vertical area centered at 0.43 MeV indicates 
the experimental centroid and width of the 7 He g.s. [HOI- m 
all calculations the lowest three 6 He states have been included 
in the binary-cluster basis. See text for further details. 

the NCSM/RGM sector of the basis. These include both 
positive-parity states, some of which are shown in Fig. [I] 
and negative-parity excitations, e.g., the 1~ soft dipole 
excitation etc. While NCSM/RGM calculations with a 
large number of clusters' excited states can become pro- 
hibitively expensive, the coupling of a few NCSM eigen- 
states of the composite system is straightforward. 

The NCSM/RGM and the NCSMC phase shifts for 
the n + 6 He five P-wave and 2 S , 1 / 2 channels are shown in 
Fig. [3] All curves have been obtained including the low- 
est three 6 He states. The NCSMC calculations (panel 
(b)) additionally include ten 7 He NCSM eigenstates, as 
described above. As expected from a variational calcu- 
lation, the introduction of the additional A-body basis 
states lAXJ^T) lowers the centroid values for all 7 He 
resonances when going from NCSM/RGM (panel (a)) to 
NCSMC (panel (b)). In particular, the 7 He 3/2" g.s. and 
5/2~ excited state are pushed toward the 6 He+rt thresh- 
old, closer to their respective experimental positions. 

The experimental accepted values for the resonance 
centroids in 7 He and the possible l/2~ states are 
shown in Table IIII1 together with our calculations. For 
NCSM/RGM and NCSMC, the resonance centroids E R 
are obtained as the values of the kinetic energy in the 
center of mass for which the first derivative of the phase 
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experiment 




NCSMC 


NCSM/RGM 


NCSM 




e r r 


Ref. 


Er r 


Er r 


Er 


3/2" 


0.430(3) 0.182(5) 


® 


0.71 0.30 


1.39 0.46 


1.30 


5/2" 


3.35(10) 1.99(17) 


[40] 


3.13 1.07 


4.00 1.75 


4.56 


1/2- 


3.03(10) 2 


M 


2.39 2.89 


2.66 3.02 


3.26 




3.53 10 


[15] 










1.0(1) 0.75(8) 


[5] 









TABLE III: Experimental and theoretical resonance centroids 
and widths in MeV for the 3/2 _ g.s. , 5/2 _ and 1/2 - excited 
states of 7 He. See the text for more details. 



shifts is maximal [4l|. The resonance widths are then 
computed from the phase shifts according to (see, e.g., 
Ref. E2) 



r = 



A8{E km )/&E k 



(4) 



An alternative, less general, choice for the resonance en- 
ergy Er could be the kinetic energy corresponding to a 
phase shift of ir/2 (thin dashed lines in Fig. [3]). While 
Eq. ((U) is safely applicable to sharp resonances, broad 
resonances would require an analysis of the scattering 
matrix in the complex plane. As we are more interested 
in a qualitative discussion of the results, we use here the 
above extraction procedure for broad resonances as well. 
The two alternative ways of choosing Er lead to basi- 
cally identical results for the calculated 3/2]" resonances, 
however the same is not true for the broader 5/2~ and 
the very broad 1/2 - resonances. The it/ 2 condition, par- 
ticularly questionable for broad resonances, would result 
in E R ~ 3.7 MeV and T ~ 2.4 MeV for the 5/2" and 
E R ~ 4 MeV (see Fig. O and T - 13 MeV for the 1/2" 
resonance, respectively. 

The resonance position and width of our NCSMC 3/2 _ 
g.s. slightly overestimate the measurements, whereas the 
prediction for the 5/2 _ is lower compared to experi- 
ment 0, [4(|, although our determination of the width 
should be taken with some caution in this case. As for 
the l/2~ resonance, the experimental situation is not 
clear as discussed in the introduction and documented 
in Table [TTT1 While the centroid energies of Refs. 11, 12] 



and 15[ are comparable, the widths are very different. 
With our determination of Er and V, the NCSMC re- 
sults are in fair agreement with the neutron pick-up and 
proton-removal reactions experiments ll|, Il2j | and defi- 
nitely do not support the hypothesis of a low lying (Er~1 
MeV) narrow (r < 1 MeV) 1/2 - resonance [4|48j]. In ad- 
dition, our NCSMC calculations predict two broad 6 P^/2 
resonances (from the coupling to the two respective 6 He 
2+ states) at about 3.7 MeV and 6.5 MeV with widths of 
2.8 and 4.3 MeV, respectively. The corresponding eigen- 
phase shifts do not reach 7r/2, see Fig. [31 In experiment, 



there is a resonance of undetermined spin and parity at 
6.2(3) MeV with a width of 4(1) MeV 0. Finally, it 
should be noted that our calculated NCSMC ground state 
resonance energy, 0.71 MeV, is lower but still compatible 
with the extrapolated NCSM value of 0.98(29) MeV (see 
Tables U and |TTT|. 

In conclusion, we introduced a new unified approach to 
nuclear bound and continuum states based on the cou- 
pling of the no-core shell model with the no-core shell 
model/resonating group method. We demonstrated the 
potential of the NCSMC in calculations of 7 He reso- 
nances. Our calculations do not support the hypothesis 
of a low lying 1/2 - resonance in 7 He. 
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